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Abstract 

We discuss gauge coupling unification in models with additional 1 to 4 complete vector-like 
families, and derive simple rules for masses of vector-like fermions required for exact gauge coupling 
unification. These mass rules and the classification scheme are generalized to an arbitrary extension 
of the standard model. We focus on scenarios with 3 or more vector-like families in which the values 
of gauge couplings at the electroweak scale are highly insensitive to the grand unification scale, 
the unified gauge coupling, and the masses of vector-like fermions. Their observed values can be 
mostly understood from infrared fixed point behavior. With respect to sensitivity to fundamental 
parameters, the model with 3 extra vector-like families stands out. It requires vector-like fermions 
with masses of order 1 TeV - 100 TeV, and thus at least part of the spectrum may be within the 
reach of the LHC. The constraints on proton lifetime can be easily satisfied in these models since 
the best motivated grand unification scale is at ~10^^ GeV. The Higgs quartic coupling remains 
positive all the way to the grand unification scale, and thus the electroweak minimum of the Higgs 
potential is stable. 



I. INTRODUCTION 



Models for new physics at the TeV scale are typically motivated by the hierarchy problem. 
They strive to explain the hierarchy between the electroweak (EW) scale and the Planck 
scale, or at least remove the incredible fine tuning required in the standard model (SM) for 
having such a hierarchy. However, the SM is stubbornly surviving the first tests at the LHC 
and there are no traces of new physics yet. In addition, the mass of the recently discovered 
Higgs-like particle suggests that the SM can be a consistent theory all the way to the Planck 
scale. This gives more weight to speculations that there is no mechanism, no new physics, 
that stabilizes the hierarchy, or that the EW scale is selected based on anthropic reasoning. 

However, even when we ignore the hierarchy problem, the SM is still not very satisfactory. 
The three gauge couplings, all couplings of the Higgs boson to fermions, the Higgs mass, 
and the Higgs quartic coupling are free parameters. This motivates us to explore extensions 
of the standard model in which at least some of these parameters could be understood. 

We have recently showed that extending the standard model by three complete vector- 
like families (SM+3VFs) with masses of order 1 TeV - 100 TeV allows for the unification of 
gauge couplings [1]. Predictions for gauge couplings at the EW scale are highly insensitive 
to fundamental parameters: the grand unification scale, the unified gauge coupling, and 
the masses of vector-like fermions. Their observed values can be mostly understood from 
infrared fixed point behavior. 

In this paper we discuss gauge coupling unification in detail in models with additional 
1 to 4 complete vector-like families (VPs), and derive simple rules for masses of vector-like 
fermions required for exact gauge coupling unification. We then focus on scenarios with 3 or 
more vector-like families that lead to insensitive unification of gauge couplings. Requiring 
the smallest splitting between masses of vector-like fermions we show that the best motivated 
grand unified theory (GUT) scale is at ~10^^ GeV. We provide examples of the spectrum 
as a function of the GUT scale, which can be as large as the Planck scale. We discuss 
constraints from proton decay and show that predictions from the best motivated region are 
close to current limits. However, due to insensitivity of predicted EW scale values of gauge 
couplings to GUT scale parameters, no sharp predictions can be made without knowing the 
spectrum of vector-like fermions. 

The focus on complete families follows from the fact that quantum numbers of quarks and 
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leptons in the SM nicely fill representations of a GUT symmetry, 10 and 5 of SU{5) or 16 
of 5*0(10). This provides a support for the idea of grand unification and the unification of 
gauge couplings [2]. Additional complete families represent some of the simplest extensions 
of the SM that can be embedded into simple GUTs.^ Consequently, there are many studies 
exploring various features of vector-like families (mostly in supersymmetric models), see for 
example Refs. [IHH]. 

In addition, vector-like fermions, not necessarily coming in complete GUT multiplets, 
are often introduced on purely phenomenological grounds to explain various discrepancies 
between observations and SM predictions. Examples include discrepancies in precision EW 
Z-pole observables p4T2]. and the muon g-2 anomaly [13]. However, with arbitrary new 
particles there are many possibilities for gauge coupling unification.^ Therefore, we gen- 
eralize the mass rules and the method to classify scenarios consistent with gauge coupling 
unification to an arbitrary extension of the standard model. 

The method to classify scenarios consistent with gauge coupling unification in terms of 
physical masses of extra particles starts with finding the mass scales that represent "average" 
masses of all particles charged under given gauge symmetry required for gauge coupling uni- 
fication (they are defined precisely in the next section and are referred to as crossing scales). 
These crossing scales are easy to obtain and they immediately give us information about 
the required spectrum. First of all, if they do not exist between the EW scale and the GUT 
scale, the gauge coupling unification in a given model is not possible, no matter what the 
splitting between masses of extra particles is. Second, the splitting between crossing scales 
represents the minimum necessary splitting in the spectrum required. Third, from the mass 
formulas that define crossing scales in terms of masses of extra particles one can immedi- 
ately see the basic features of the spectrum required, and the spectrum can be calculated. 
In addition, these formulas also indicate the freedom one has in imposing further relations 
between masses of extra particles. This might be useful when searching for models that 

^ This does not mean that the masses of vector-like fermions needed for gauge coupUng unification necessar- 
ily result from a simple unified boundary condition. By simple GUTs we mean that there is no additional 
mechanism required to keep particles in incomplete GUT multiplets significantly below the GUT scale, or 
to split their masses over many orders of magnitudes that would, to large extent, ameliorate the motivation 
for GUTs. 

^ For examples of recent studies investigating the effects of extra particles on gauge coupling unification in 
models without supersymmetry, see Refs. |14[ I15j. 



3 



relate masses of particles at a given scale. The mass rules given in terms of particle masses 
can be evolved to an arbitrary scale, e.g. the GUT scale, which would provide the boundary 
conditions that need to be satisfied. However, the renormalization group (RG) evolution of 
the mass rules depends on additional assumptions one has to make about the origin of the 
masses and the scale at which these masses are generated.^ 

This paper is organized as follows. In Sec. |TT]we discuss RG evolution of gauge couplings in 
models with extra VFs. We start with the discussion of IR fixed point predictions for gauge 
couplings, then add threshold corrections from a universal mass of vector-like fermions, and, 
finally, we add effects from splitting masses of vector-like fermions. We discuss sensitivity 
of predicted values of gage couplings to fundamental parameters. Finally, we derive simple 
mass rules that have to be satisfied in order to get exact gauge coupling unification. We 
generalize the method to classify all solutions consistent with gauge coupling unification to 
an arbitrary extension of the SM. In Sec. |III| we discuss constraints from proton decay, the 
stability of the EW minimum of the Higgs potential, and discuss a possible origin of masses 
of vector-like fermions. We give few concluding remarks in Sec. |IV[ 

II. RENORMALIZATION GROUP EVOLUTION OF GAUGE COUPLINGS 

The one-loop renormalization group equations (RGBs) for three gauge couplings, = 
gf/4:'7i, are given by: 

^ =/?(«,) = (1) 

dt ^ ^ 27r ' ^ ^ 

where t = IylQ/Qq with Q representing the energy scale at which gauge couplings are 
evaluated. The beta function coefficients, bi, in the SM with nf families are given by 

/ 1 4 43 4 4 \ 

bi = h -nf, h -Uf, -11 + -Uf . (2) 

VIO 3 ^' 6 3 ^' 3 V 

For nf = 3 we get the usual SM result, 6j = (41/10, —19/6, —7). With extra pairs of 
complete VFs we have nj = 3 + 2xA^(a vector-like partner contributes in the same way). 

^ The study of gauge coupling unification is, to large extent, unaffected by these assumptions; only the 
physical masses of particles matter in the leading order. If the masses originate from Yukawa couplings 
to extra scalars that get vacuum expectation values at an intermediate scale, these may contribute to the 
RG evolution of gauge couplings at 2-loop level. However unless the extra couplings are large these effects 
would be negligible. 
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FIG. 1: RG evolution of gauge couplings: 03 (top solid line), 02 (middle solid line), and ai (bottom 
solid line) in the SM extended by three vector-like families for ac = 0.3 at Mq = 2 x 10^® GeV. 
Dashed lines in the same order show the running of gauge couplings in the SM. Masses of 3 VFs are 
neglected in the left plot, and fixed to 10 TeV (indicated by Myp) in the right plot. The crossing 
points in the evolution of gauge couplings in the SM+3VFs and the SM indicated in the left plot 
define the common threshold scales, Mi^2,3, for masses of particles charged under given symmetry 
required for exact gauge coupling unification. 

For example, in SM+3VFs we find bi = (121/10, 29/6, +1) which indicates that all three 
gauge couplings are asymptotically divergent (this result obviously holds for 3 or more pairs 
of VFs). 

The evolution of gauge couplings in the SM and an example of the evolution in the 
SM+3 VFs case are shown in Fig. [l] The numerical analysis closely follows that of Ref . [1] . 
For the SM evolution we use the Z-scale central values of a^l^{Mz) = 127.916, sin^ Ow = 
0.2313, and a^lMz) = 0.1184, together with the top quark mass rrit = 173.2 GeV, which 
can be found in Ref. [3j. The apM and sin^ 9w are related to ai^2{Mz) through 

a' 



sin^ e 



w 



and a EM = a2sin^6'vi/, 



(3) 



a2 + a' 

where, assuming the SU{5) normalization of the hypercharge, a' = (3/5)ai. We set the 
Higgs boson mass to rrih = 126 GeV [151 [IZ] . The example of the RG evolution of gauge 
couplings in the SM+3VFs starts with unified gauge coupling ac = 0.3 at Mq = 2 x 10^^ 
GeV. The crossing points in the evolutions of gauge couplings in these two cases, which will 
be important for the discussion of threshold corrections, are indicated in the left plot by 
Mi_2,3- In all numerical results we use full two loop RGEs [18], we integrate out all particles 
with masses above Mz at their mass scale, and include one-loop matching corrections for 



rrit and ruh fL9\. We assume that Yukawa couplings of vector-like fermions are negligible, 
and we also neglect Yukawa couplings of all fermions in the SM except the top quark. 

The results of the numerical analysis we present can be understood from approximate 
analytic formulas. The one-loop RGEs can be solved, and we can express gauge couplings 
at the EW scale in terms of the GUT scale, and values of gauge couplings at Mq: 

c^7\Mz) = |^ln§^ + ai\MG). (4) 
zvr Mz 

Assuming gauge coupling unification, ai^Mc) = ac, and neglecting threshold corrections 
both at the EW scale and the GUT scale, we can express one gauge coupling in terms of 
the other two. For example: 



= (6, -fr,)4, + 3/1(63 

where = sin^ 9w{Mz), and c^r = cos'^ 9w{Mz). For the measured values of uem and 
the SU(5) embedding of the SM predicts a^i^Mz) — 0.07 which is about 40% below the 
experimental value. 

Adding complete chiral or vector-like families at the EW scale does not change at all the 
one-loop prediction given in Eq. ^ since complete families contribute equally to all three 
beta function coefficients, see Eq. ([2]). Furthermore, the scale of unification (more precisely 
the scales where any two couplings meet) does not change at one-loop, only the value of 
the unified gauge coupling increases. With increasing the number of extra families, at some 
point, the couplings become non-perturbative before they meet, and eventually reach the 
Landau pole. Further increase of the number of families lowers the energy scale at which 
the Landau pole occurs. 

However, the SM extended with a sufficient number of complete vector-like families so 
that all couplings are asymptotically divergent offers a new possibility. Vector-like families 
introduce an additional scale to the problem associated with masses of vector-like fermions, 
MvF, and they contribute to the RG evolution of gauge couplings only above this energy 
scale. This allows us to consider models with a large (but still perturbative) unified gauge 
coupling at a high scale, higher than the scale at which the Landau pole would occur if 
the VFs were at the EW scale. Consequently, in the RG evolution to lower energies, gauge 
couplings run to the (trivial) infrared (IR) fixed point. Thus, at lower energies, the values of 
gauge couplings are determined only by the particle content of the theory and how far from 



the GUT scale we measure them. Since the exact value of ac becomes irrelevant, instead of 
one prediction of the conventional unification, Eq. (|5]), we have two predictions for ratios of 
gauge couplings. At the Myp scale, the vector-like fermions are integrated out, and below 
this scale gauge couplings run according to the usual RG equations of the standard model. 
In a way, the two parameters of the conventional unification, Mq and ac-, are replaced by 
Mg and Mvf- The discrepancies of IR fixed point predictions from observed values can be 
explained by threshold effects of extra vector-like fermions. 



A. IR fixed point predictions for gauge couplings 

The IR fixed point predictions were discussed in detail in Ref . p] . In models with asymp- 
totically divergent couplings, these can be easily obtained if the 1-loop RGEs are good 
approximations. Assuming a large enough unification scale and large (but still perturba- 
tive) unified gauge coupling, the first term in Eq. (|4]) dominates, and the ratios of gauge 
couplings are given by ratios of beta function coefficients. 



(6) 



a^{Mz) _ 6j 
aj{Mz) ~ hi' 

This can be translated into the prediction for sin^ 6w'- 

sm Ow = ^ — - = , , (7) 

where h' = (5/3)6i. Numerically, we find sin^ 9w = 0.193 in the case of SM+3VFs, which is 
identical to the value obtained assuming 9 chiral families [201 EI]- Similarly, in SM+4VFs 
we find sin^ Ow = 0.234. 

In the case of SM+3VFs, the one-loop RGE for as given in Eq. ([T]) is not a good approx- 
imation because of the accidentally small 63 coefficient. The two-loop contribution to the 
beta function is well approximated by the term proportional to 

dc(3 ^/ N , «3 /^x 

(8) 

where B3 = -102 + (76/3)n/ = 126 for SM+3VFs [E]. Thus, the two loop contribution is 
larger than the one- loop contribution for ^3 > 0.1.^ 

^ This is a consequence of a very small 1-loop beta function coefficient and it is not an indication of non- 
perturbativity. The coupling is still perturbative, the dominant 3-loop contribution to /^(aa), proportional 
to Q?!, represents a ~ 5% correction to 1 + 2-loop beta function for 03 ~ 0.1 [52]. 
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The RGE for 0:3 can be solved by adding the 1-loop contribution as an expansion in e = 
Anb^i/B^ to the solution obtained from the 2-loop contribution only [23], PQ. Alternatively, 
we can solve the full RGE given in Eq. ([8| and find: 



a^\Mz) --ln(l + —^] = ^\n^ + a7,'--\n(l + ^]. (9) 



Neglecting a^^, we obtain the second prediction: 

az{Mz) _ h2 + h' 



aEM{Mz). (10) 



Numerically, for aEhii^Mz) = 1/127.916, it predicts a^i^Mz) — 0.072 in the case of 
SM+SVFs. 

The beta function coefficients for 0:3 in the SM+4VFs scenario are 63 = 11/3 and B3 = 
530/3. The 1-loop term in the RG equation ^ dominates for a-^ < 0.26 in this case. 

The proximity of predictions from the IR fixed point, Eqs. ([T]) and (10), to observed 
values is certainly intriguing. Although they are not a perfect match to measured values, 
the discrepancies can be easily accommodated by taking into account threshold corrections 
from vector-like fermions that should be integrated out at the Mvf scale. 



B. Mass scale of vector-like fermions and sensitivity to fundamental parameters 

The existence of a scale associated with masses of vector-like fermions is strongly sug- 
gested by the overlay of the RG evolution of gauge couplings in the SM and those in the 
SM+3VFs assuming unified gauge coupling at a high scale given in Fig. [TJ All three gauge 
couplings in these two scenarios cross at comparable scales suggesting a common threshold 
at which particles from VFs are integrated out. Indeed, for the example given in Fig. [T| 
fixing all the masses of 3 VFs to 10 TeV, shown in Fig. [l] (right), the EW scale values of 
gauge couplings are predicted within 8% from measured values. In the next subsection, 
we will show that the measured values of gauge couplings can be precisely reproduced by 
splitting the masses of vector-like fermions. First, however, we would like to discuss general 
features of this result assuming the common mass of VFs. 

The fairly good agreement of predicted values of gauge couplings from 3 VFs at ~ 10 
TeV with observed values does not rely on the specific choice of the GUT scale and the value 
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of the unified gauge coupling. The EW scale values of gauge couplings are highly insensitive 
to these parameters, which can be understood from IR fixed point behavior. 

The low sensitivity of predicted values of gauge couplings to fundamental parameters 
is demonstrated in Fig. [2] (left). It shows a large region of the GUT scale, Mg, and the 
universal mass of fermions from 3 vector-like families. My p., from which the values of gauge 
couplings at the EW scale are simultaneously predicted within 10% from the measured 
values. It also shows the best fit that predicts all couplings within 6%. The GUT scale is 
the best motivated between 10^^ GeV and 10^'^ GeV with the best fit close to 10^^ GeV. For 
completeness, a similar plot for aEMt -s^, and 0:3 is presented in Fig. |2] (right). However, as 
we will see from the discussion of threshold corrections in the next subsection, the plot on 
the left for ai, 02, and 03 is more indicative of the best motivated values of Mq and My p. 

In order to understand the sensitivity of the EW scale values of gauge couplings to 
fundamental parameters quantitatively, it is instructive to estimate separate contributions 
to ai,2,3(^z) from ac M^, and My p. Note that the values of 0^23(^2) are approximately 
59, 30, and 8.4, respectively. From Eqs. ^ and (10), we see that ac ^ 0.3 contributes less 



than ~ 10% to the EW scale values of gauge couplings. It is the least important parameter. 
Plots in Fig. |2] for any ac > 0.3 would look almost identical. Increasing ac moves all 
the contours slightly to the right. The largest contribution to EW scale values of couplings 
originates from the 1 / {2tt) \ii{Mc/ Alz) — 5.2 term multiplied by corresponding beta function 
coefficients. 

The second largest contribution to gauge couplings comes from masses of vector-like 
fermions. The IR fixed point predictions for the gauge couplings at the EW scale, obtained 
from Eq. ^ for ai 2 with ai 2(Mg) = ac, and from Eq. ^ for 03, are modified by threshold 
corrections Tj: 

1 - ai{Mz)Ii 

that depend on masses of the extra vector-like fermions. These threshold effects are well 
approximated by the leading logarithmic corrections:^ 

T, = -yb{\n^, (12) 

f 

^ These corrections correspond to removing one loop contributions of vector-like fermions from Eqs. Q and 
Eq. ([9| below their mass. It is an excellent approximation for ai.2 and sufficient approximation for 
since, for the IR value of as, the 1-loop term in the RG equation dominates. 
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FIG. 2: Left: contours of constant values of predicted gauge couplings at Mz, ai (green), a2 
(blue), and 03 (red), as functions of the GUT scale, Mq, and the universal mass of fermions from 
3 vector- like families, Myp, for fixed ac = 0.3. Solid lines represent the central experimental 
values of three gauge couplings, the shaded regions represent ±10% ranges, and the dashed lines in 
unshaded areas represent ±20% ranges. The lightly shaded area corresponds to a ±50% range of 
03. In the overlapping (bright red) region, all three gauge couplings are simultaneously predicted 
within 10% from the measured values, and the small black area in the red region represents the 
best fit with all three couplings within 6% from the measured values. The gray region corresponds 
to a2,{Mz) > 0.3; a3(M^) becomes non-perturbative very fast with increasing Myp from the value 
that corresponds to the boundary of this region. Right: the same as in the plot on the left but for 
0(EM (orange), (purple), and 03 (red). The dotted purple line represents s^r being —5% from 
the central value. 

where b{ is the contribution of a given fermion /, with mass Mf, to the corresponding beta 
function coefficient [IB]. For particles originating from vector-like families, these contribu- 
tions, summarized in Table |T[ are identical to contributions from fermions in the standard 
model. The contribution from the complete family is identical to all three beta function 
coefficients and equal to 4/3 for a chiral family, and 8/3 for a vector-like pair (16 ± 16 in 
the SO(IO) language). 

The correction to of about ±40% is crucial in order to reproduce the measured value. 
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TABLE I: Quantum numbers and contributions to beta function coefficients of particles from extra 
vector-like families. The names are chosen to mimic those of the standard model particles with 
the same quantum numbers. For each particle, there is a corresponding vector-like partner, and 
its contributions to the beta function coefficients are identical. The bi coefficients correspond to 
the SU(5) normalization of the hypercharge. 



Particle 


SU{3) X 


SU{2) 


X C/(l) 


^3 


b2 


bi 


Q 


3 


2 


1/6 


2/3 


1 


1/15 


u 


3 


1 


-2/3 


1/3 





8/15 


E 


1 


1 


1 








2/5 


L 


1 


2 


-1/2 





1/3 


1/5 


D 


3 


1 


1/3 


1/3 





2/15 



As can be seen in Figs, m and |2| it is indeed that determines Myp — 10^ GeV and, 
consequently, Mg — 10^^ GeV. The other two couphngs are within 10% from measured 
values in much larger ranges of Mvf and would actually prefer smaller Mvf and Mg- 

Out of the three parameters, the EW scale values of gauge couplings are the most sensitive 
to changes in My p. However, since Myp is only responsible for at most ~ 40% of the EW 
scale values of couplings, the overall sensitivity is still very small. Most of the EW scale 
values of couplings originate from the IR fixed point. Since no precise cancellations between 
separate contributions are required, there are large ranges of fundamental parameters from 
which the predicted values of gauge couplings at the EW scale are close to observed values. 

The standard model extended by 4 vector-like families (SM+4VFs) allows for insensitive 
unification of gauge couplings in a similar way as the SM+3VF. Predicted values of gauge 
couplings at Mz as functions of the GUT scale and the universal mass of fermions from 
4 vector-like families for fixed og = 0.3 are shown in Fig. [3j There are, however, notable 
differences from the SM+3VFs case. First of all, the common mass of vector-like families 
moves to ~ 10^ — 10'^ GeV. This is easily understood from the fact that more matter makes 
gauge couplings run faster, and thus the VFs must stop contributing to RG evolution at a 
higher scale; otherwise, the EW scale values of gauge couplings would be too small. Second 
of all, the 1-loop IR fixed point value of sin^ 9iy is 0.234 which is larger than in the SM+3VFs 
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FIG. 3: The same as in Fig. [2] but for the SM extended by 4 vector-Uke famihes. 

case and actually very close to the measured value. Overall, this however does not make 
the predictions much better than in the SM+3VFs case, since requires Myp larger than 
the one needed to reach the measured value of sin^ 6^. Finally, as a result of larger masses 
of VFs required in the SM+4VFs scenario, the sensitivity of EW scale values of gauge 
couplings to fundamental parameters increased, which is visible in Fig. |3] as narrower 10% 
bands compared to those in Fig. [2] corresponding to the case of SM+3VFs. 

It is easy to extrapolate to a larger number of VFs. Increasing the number of VFs requires 
larger Myp closer and closer to the GUT scale. The sensitivity of predicted values of gauge 
couplings to fundamental parameters is increasing and approaching the sensitivity in the 
SM. 

In the SM extended by 1 or 2 vector-like families the predictive power is lost, since the 
unified gauge coupling is small and its specific value is crucial for predictions for gauge 
couplings at the EW scale in a similar way as in the SM. The difference form the SM is that 
the exact unification of gauge couplings is now possible with split masses of VFs. We will 
include these solutions as a curiosity in the next subsection. 
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C. Threshold effects of vector-like fermions 



Let us now turn our attention to precise predictions for gauge couplings rather than a 
~10% agreement. For this we need to consider threshold effects from splitting masses of 
VFs. 

The necessity to split masses of particles from extra 3VFs is indicated in Fig. [T] (left) 
by slightly different scales at which the RG evolutions of gauge couplings in the SM and 
SM+3VFs cross. For the example in this figure the crossing scales for ai, and are 
Ml ~ 100 TeV, M2 — 1 TeV, and M3 ^ 10 TeV. These scales determine threshold corrections 



Tj = (4/7r) hi{Mi/ Mz)-, see Eq. (12), required for gauge coupling unification. Any spectrum 
that leads to required threshold corrections will reproduce the measured values of gauge 
couplings. 

The crossing scales are increasing with increasing Mq and depend very little on ac for 
ac ^ 0.3. For different values of M^, they can be read out of Fig. |2] (left) as corresponding 
values of Mvf for which we obtain the measured value of given gauge coupling. Similarly, 
in the case of SM+4VFs, the values of Mi 2,3 can be read out of Fig. |3] (left). For values of 
Mg not shown, or for other scenarios, the crossing scales can be easily calculated from RG 
equations as functions of Mq and ac- 

In general, for pairs of vector-like families, once we know values of crossing scales Mi_2,3 
for chosen GUT scale, the masses of fermions must satisfy: 

N 



AN, M3 l>r^/,o. Mo- ,rr, Mu- Md-\ 



In 



■t=i 

N 



4iV, M2 lv^/,o, Mq. Ml.\ 

i--^ = -E^?i-i;^ + ^^i-i;^ ' (14) 



37r Mz TT ^ V Mz Mz 

i=l 

4iV,__ Ml 1 ^ Mq^ , Mu, , Md, , Mr , Mj, 



In 



TT ^ V Mz Mz Mz Mz Mz J 



37r Mz vr 

i=l 

in order to get exact gauge coupling unification at given GUT scale. In the case of universal 
masses of particles with the same quantum numbers, e.g. Mq^ = Mq^ = ■ ■ ■ = Mqj^ = Mq, 
these mass rules can be written in a simple form: 

m//' = n ^=1,2,3. (16) 

F=Q,U,D,L,E 



13 



Inserting the beta function coefficients from Table |H we find: 

M| = M^MuMd, (17) 
= M^Ml, (18) 
= MQMfjMlMlM%. (19) 

These formulas hold for any number of complete vector-like families, only the values of Mi 2,3 
depend on the specific scenario. In the case of non-universal masses of particles with the 
same quantum numbers the above formulas are still valid with the replacement: 

Mf = [Mf.Mf, . . . Mp^f , F = Q,U,D, L, E. (20) 



From Eqs. (17) - (20) we can immediately see that splitting fermions with the same 
quantum numbers does not help to find a solution if a solution does not exist with universal 
masses. Thus, it is sufficient to assume universal masses, Mp, of particles with the same 



quantum numbers, and Eqs. (17) - (19) classify possible solutions. In addition, for each 
solution with universal masses, there are other solutions with split masses, and the only 
constraint is that their geometric mean is the universal mass needed for a given solution. 
There are many solutions available, since we have 5 different masses that have to satisfy 



three conditions (17) - (19) . However, it is not guaranteed that for a given GUT scale there 
is a phenomenologically viable solution. Clearly, the crossing scales have to be above the 
EW scale, and even then the solution might require new fermions below experimental limits 
or some fermions above the GUT scale or the Planck scale. 

Representative examples of the spectrum for various values of Mq in the case of SM+3VFs 
are given in Fig. |4j The value of ac is fixed to 0.3; however, the spectrum is not very sensitive 
to this choice as previously discussed. The spectrum shown is just an example, motivated by 
the smallest splitting between masses required, and it is not unique. A specific example with 
exact numerical values was also given in Ref. [1]. The GUT scale motivated by the lowest 
splitting required between masses of vector-like fermions is at ~ 10^^ GeV in agreement with 
what is suggested in Fig. [2] (left), and the masses are split between ~1 TeV and ~100 TeV. 

There is a lower bound on the possible GUT scale at ~ 10^^ GeV. For smaller Mq the 
crossing scale M2 is too small, see Fig. [2] (left), and thus a phenomenologically viable solution 
does not exist. With increasing Mg, the splitting of fermion masses is increasing, which can 
also be inferred from a larger splitting of crossing scales. The GUT scale can be as high as 
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FIG. 4: Masses of vector-like fermions leading to exact gauge coupling unification as functions of 
the GUT scale in the case of SM+SVFs. The universal mass for particles with the same quantum 
numbers is assumed. The value of ac is fixed to 0.3. Smaller values of Mq (not shown) can still be 
consistent with gauge coupling unification for smaller ac- The spectrum shown is just an example, 
it is not unique. 

the Planck scale. However, in that case the masses of vector-like fermions have to be split 
over 6 orders of magnitude. 

Note that quark doublets, Q, are typically predicted at ~1 TeV. The preference for Q 
being the lightest of vector-like fermions can be understood from M2 < Mi 3. However, 
there are also solutions with L being the lightest. Keep in mind, however, that these 
masses represent geometric means of masses of particles with the same quantum numbers. 
Therefore, when considering split masses of fermions with the same quantum numbers, any 
fermion can be the lightest one and as light as current experimental limits. 

For SM-|-4VFs, examples of the spectrum are given in Fig. |5} The main features are 
very similar to the case of SM+3VFs. The GUT scale motivated by the lowest splitting 
required between masses of vector-like fermions is also at ~ 10^^ GeV, in agreement with 
what is suggested in Fig. [s] (left), and about two orders of magnitude splitting of masses of 
vector-like fermions is required. The main difference from the SM+3VFs case is that the 
spectrum shifted to 10^ - 10^ GeV. 

For completeness, we also include examples of the spectrum needed for exact gauge 
coupling unification in the case of SM+IVF in Fig. [6] and SM-|-2VFs in Fig. [7} For these 
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FIG. 5: The same as in Fig. [4] but in the case of SM+4VFs. 
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FIG. 6: The same as in Fig.|4]but in the case of SM+IVF. In this case, values of ac are optimized 
for given GUT scale, and are close to 0.03 for all Mq shown. For smaller or larger values of Mq, 
the unification is not possible for any spectrum. 

cases, the EW scale values of gauge couplings are highly sensitive to ac- Thus ac in these 
examples is not fixed but rather optimized for the given Mq- In both cases, the exact 
unification can be achieved even in the region consistent with limits on proton lifetime. 
However, the required splitting between masses of vector-like fermions is sizable. 
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FIG. 7: The same as in Fig. [4] but in the case of SM+2VF. In this case values of ac are optimized 
for given GUT scale, and vary between 0.042 and 0.048. Smaller values of Mq (not shown) can 
still be consistent with gauge coupling unification. 

D. Generalization of mass rules to other extensions of the SM 

The mass rules we have just derived can be generalized to any extension of the SM. The 
existence of crossing scales is a necessary condition for achieving gauge coupling unification 
in a given model. This follows from the fact that integrating out extra fields above the EW 
scale can only increase gauge couplings at the EW scale. Therefore, values of predicted 
couplings at the EW scale without considering the mass effect of extra matter fields have 
to be smaller than the measured values. The crossing scales depend only on ac and Mq- 
Thus, requiring that the crossing scales exist leads to limits on possible values of the GUT 
scale and ac- 

For chosen ac and Mq-, we can find the crossing scales Mi 2,3 for all three gauge couplings. 
If one loop RGEs are good approximations, these crossing scales can be easily found by ap- 
plying Eq. ^ separately between Mz the Mj scales using the SM beta function coefficients, 
6f^, starting with the observed values of Q.i^exp{,Mz), and between Mi and Mq scales using 
beta function coefficients in the given extension, hi, assuming gauge couplings exactly unify. 
We get: 

Mi 27r f . _i hi , Mg 



= [^o-i.^^ + + ^ 1" ^ j . (21) 

The meaning of crossing scales is the same as in extensions of the SM with VFs; namely, 
they represent the threshold corrections, Tj = {pi — hf^'^)/{2Ti)\n.{Mi/Mz), that masses of 
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extra particles must generate in order to reproduce the measured values of gauge couplings 
starting from the given ac and Mq- The rest follows what we did for complete VFs. Once 
we know values of crossing scales Mi 2,3, in order to get exact gauge coupling unification, 
the masses of extra particles must satisfy: 

Mf'"'''"^ = n^?' ^=1,2,3, (22) 

F 

where the product is over all extra fermions (or scalars) charged under a given gauge sym- 
metry. For a vector-like pair of fermions, the corresponding mass on the right-hand side 
appears twice. As in the case of complete VFs, it is sufficient to consider the universal mass 



of all particles with the same quantum numbers. The universal mass that enters Eq. (22) 
represents their geometric mean. 



For any model with an arbitrary particle content, the crossing scales (21) as functions of 



ac and Mq together with the mass rules (22) classify all the solutions consistent with gauge 
coupling unification in terms of physical masses of extra particles. 

Let us illustrate the usefulness of crossing scales and the mass rules on one example. 
Let us ask if there is any spectrum of extra particles in the SM extended by one vector- 
like family that leads to exact gauge coupling unification for Mq = 10^® GeV. This choice 
corresponds to one of the points in Fig. |6} and so we already have the answer we can compare 
with. However, this answer is obtained by a fairly complicated numerical procedure that 
iteratively solves coupled differential equations with masses of extra vector-like fermions 
varied untill the EW scale values of gauge couplings are precisely reproduced. Using our 
method, we can get the basic features of the required spectrum fast. 

For ac = 0.0286 that corresponds to the given example in Fig. [HI the crossing scales 



Mi,2,3, easily calculated from Eq. (|2l|), are 1 x lO^^ GeV, 7 x 10^ GeV, and 2 x 10*^ GeV. 
This immediately tells us that there will be more than 8 orders of magnitude splitting 
between masses required. Knowing the crossing scales, we can easily see the basic features 



of the spectrum that will work. From Eqs. (22 ), which in this case are the same as Eqs. ( 17) - 



(19), we see that Mq, which heavily weighs on M2, should be less than M2, while everything 
with large hypercharge (especially E and U) should be above Mi in order to find a solution. 
For a specific example, one can choose two masses and calculate the rest of the spectrum 



from Eqs. (17) - (19). Given a large splitting between M2 and Mi in this case, it would be 



easiest to choose the masses of Q and U as a starting point. Once we have one solution, 
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varying the starting masses of Q and U and calculating the rest of the masses from Eqs. (17) 
- (19) will give us all possible solutions for the given ac and Mq- This procedure can be 
repeated for any ac and M^, or the solutions can be plotted as functions of these variables. 



III. DISCUSSION 

So far, we have only considered constraints on the GUT scale and masses of vector-like 
fermions from gauge coupling unification. In order for this scenario to be easily embedded 
into simple grand unified theories, based on SU{5) or 5*0(10), the constraints on proton 
lifetime and the stability of the EW minimum of the Higgs potential should be satisfied. 

The most stringent limits on proton lifetime come from Super-Kamiokande. For the 
dominant decay mode from dimension-6 operators, the limit is r(p — > 7r°e+) > 1.4 x 10^"^ 
yrs [23]. Assuming naively, that the proton lifetime is Tp ~ Mq/ {a^m^^, where nip is the 
mass of the proton, this limit translates into the lower bound on the GUT scale: Mq > 
1.5 X 10^^ GeV for ac = 0.3 which we use in our examples. However, the prediction for 
the proton lifetime is somewhat model dependent (see for example Refs. [21 |2H ES] and 
references therein), and so we do not impose the strict limit in the plots we present. In 
addition, the plots would look very similar for any large value of ac, but the limits would 
differ. The interested reader can easily impose the limit on any scenario by simple rescaling 
of the mentioned limit using the formula for the proton lifetime. 

It is interesting to note that the best motivated value of the GUT scale is in the ~10^^ GeV 
range which is basically at the current limit. It is, however, not possible to make precise 
predictions without knowing the masses of vector-like fermions. For example, a scenario 
with the GUT scale larger by a factor of 3 results in ~2 orders of magnitude enhancement 
of the proton lifetime but would only require modest changes in the spectrum of vector-like 
fermions in order to have exact gauge coupling unification. This inability to make precise 
predictions of GUT scale parameters is a direct consequence of the insensitivity of the EW 
scale couplings to GUT scale boundary conditions. 

The RG evolution of the top Yukawa and Higgs quartic couplings in the SM+3VFs for 
Mg = 2 X 10^^ GeV and ac = 0.3 was given in Ref. [I]. The Higgs quartic coupling 
remains positive all the way to the GUT scale, and thus the electroweak minimum of the 
Higgs potential is stable. This result holds in a large range of Mg and ac, especially in the 
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best motivated region. Therefore, these scenarios represent some of the simplest possible 
extensions of the standard model that can be embedded into grand unified theories, with a 
sufficiently long lived proton, and a stable EW minimum of the Higgs potential. 

We have not investigated the origin of masses of vector-like fermions needed for gauge 
coupling unification. This would require additional assumptions about the mechanism that 
generates them and the scale at which boundary conditions are set. The masses of vector- 
like fermions may be fundamental lagrangian parameters, or they can originate from Yukawa 
couplings to one or several additional scalars (singlets under SM gauge symmetry, but possi- 
bly charged under family symmetries) that acquire vacuum expectation values at any scale 
between the GUT scale and the Myp scale. In addition, vector-like fermions can have non- 
zero Yukawa couplings to the SM Higgs doublet, which add another layer of complexity 
by contributing to the physical masses and possibly significantly affecting the RG evolu- 
tion of other parameters that directly determine their masses. The study of gauge coupling 
unification is, to a large extent, unaffected by these assumptions, only the physical masses 
of particles matter in the leading order. Fundamental lagrangian masses would not affect 
the running of gauge couplings at all, and the Yukawa couplings to extra scalars may only 
contribute to the RG evolution of gauge couplings at 2-loop level. However, in any specific 



scenario, the mass rules (17) - (19) can be evolved to the GUT scale (or other relevant scale), 
and the freedom to choose some of the masses can be used to search for simple boundary 
conditions that are consistent with gauge coupling unification. 

Finally, it is intriguing to consider a connection with the anthropic solution to the hi- 
erarchy problem, or the EW scale [261 EZI- Adding VFs to the SM makes this possibility 
more appealing, since in the SM, special values of gauge couplings either at the EW scale or 
some high scale have to be selected. In scenarios that we discussed, the EW scale values of 
gauge couplings close to the observed values are not very special, but rather quite a generic 
outcome from large ranges of fundamental parameters. For example, in the SM+3VFs case, 
as far as Myp < 25 TeV, for any ac ^ 0.2, and Mq anywhere between 10^^ GeV and the 
Planck scale, the predicted values of gauge coupling at the EW scale are always within 50% 
of the measured values (ai,2 typically well within 20%). This is indicated by the lightly 
shaded region in Fig. |2} and a similar region is indicated in Fig. |3] for SM+4VFs. Further- 
more, if the EW scale and Mvf have the same origin, it would also explain the proximity 
of the QCD scale to the EW scale. The beta function of changes the sign at Myp, and 
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below this scale, it starts running fast toward Aqcd- 

IV. CONCLUSIONS 

We have discussed gauge coupling unification in models with additional 1 to 4 complete 
vector-like families. In scenarios with 3 or more vector-like families the values of gauge 
couplings at the electroweak scale are highly insensitive to the grand unification scale, the 
unified gauge coupling, and the masses of vector-like fermions. Their observed values can 
be mostly understood from infrared fixed point behavior. Starting with a large (but still 
perturbative) unified gauge coupling at a high scale, the values of gauge couplings at lower 
energies are determined only by the particle content of the theory and how far from the 
GUT scale we measure them. Since the exact value of ac becomes irrelevant, instead of 
one prediction of the conventional unification, we have two predictions for ratios of gauge 
couplings. These predictions are modified at the Mvf scale, where the vector-like fermions 
are integrated out, and below this scale, gauge couplings run according to the usual RG 
equations of the standard model. 

Assuming first a common mass of vector- like fermions, Myp, we showed predictions for 
three gauge couplings at the EW scale as functions of Mq and My p. We found that the 
observed values of gauge coupling are reproduced with good precision from a large range 
of parameters. Especially, Mq can be varied over several orders of magnitude while having 
all three gauge couplings within 20% from observed values. The best fit, which predicts all 
three couphngs within 6% from measured values, suggests Mq ~ lO-*^^ GeV, Mvf — 10^ 
GeV in the case of SM+3VFs, and Mvf ^ 10^ - 10^ GeV in the case of SM+4VFs. 

The best motivated GUT scale, ~10^^ GeV, predicts a proton lifetime close to current 
limits. However, due to insensitivity of the predicted EW scale values of gauge couplings to 
GUT scale parameters, no sharp predictions can be made without knowing the spectrum of 
vector-like fermions. In addition, it was previously shown that the Higgs quartic coupling 
remains positive all the way to the GUT scale, and thus the electroweak minimum of the 
Higgs potential is stable. This result holds in a large range of Mq and ac, especially in the 
best motivated region. Therefore, these scenarios represent some of the simplest possible 
extensions of the standard model that can be embedded into grand unified theories, with a 
sufficiently long lived proton, and the stable EW minimum of the Higgs potential. 
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The discrepancies of IR fixed point predictions from observed values can be explained by 
threshold effects of extra vector-like fermions. We showed examples of the spectrum for the 
GUT scale varied between 10^^ GeV and 10^^ GeV. We derived simple rules for masses of 
vector-like fermions required for exact gauge coupling unification. In addition, we generalized 
the mass rules and the method of using crossing scales of evolutions of gauge couplings in the 
SM and the given extension to classify scenarios consistent with gauge coupling unification 
to an arbitrary extension of the standard model. The problem of finding all possible mass 
spectra in a given model consistent with gauge coupling unification is reduced to solving a 
set of simple algebraic equations that masses of extra particles have to satisfy. 

With respect to the sensitivity to fundamental parameters, the model with 3 extra vector- 
like families stands out. In the best motivated region, it requires vector-like fermions with 
masses of order 1 TeV - 100 TeV, and thus at least part of the spectrum may be within the 
reach of the LHC. Notably, quark doublets, Q, are typically predicted at ~1 TeV. However, 
only geometric means of masses of particles with the same quantum numbers are constrained 
by gauge coupling unification. Therefore, when considering split masses of fermions with 
the same quantum numbers, any fermion can be the lightest one and as light as current 
experimental limits. Besides direct production of these particles at the LHC, it may be 
also possible to observe their effects in a variety of processes. However, they typically affect 
standard model predictions only through mixing with light fermions, which is highly model 
dependent. The discussion of gauge coupling unification that we focused on here is neghgibly 
affected by such mixing. 
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